Theory of carrier transport in bilayer graphene 
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We develop a theory for density, disorder, and temperature dependent electrical conductivity of 
bilayer graphene in the presence of long-range charged impurity scattering as well as an additional 
short-range disorder of independent origin, establishing that both scattering mechanisms contribute 
significantly to determining bilayer transport properties. We find that although strong screening 
properties of bilayer graphene lead to qualitative difi'erences with the corresponding single layer 
situation, both systems exhibit the linearly density dependent conductivity at high density and the 
minimum graphene conductivity behavior around the charge neutrahty point due to the formation 
of inhomogeneous electron-hole puddles induced by the random charged impurity centers. 

PACS numbers: 81.05.Uw; 72.10.-d, 73.40.-c, 72.20.Dp 



Ever since the discovery of graphene [ij its transport 
properties as functions of carrier density and tempera- 
ture have been of key fundamental and technological in- 
terest. The fundamental interest arises from the unique 
Hnear massless chiral Dirac dispersion of electrons (holes) 
in the graphene conduction (valence) band with the sys- 
tem being a zero-gap semiconductor which on doping (or 
external field induced gating) changes its character con- 
tinuously from being an 'electron-metal' to a 'hole-metal' 
as it goes through the charge neutral Dirac point. Unique 
transport properties of massless, gapless, and chiral Dirac 
particles in 2D single layer graphene (SLG) as functions 
of their density and temperature have attracted a great 
deal of experimental and theoretical attention over the 
last few years. 

More recently, however, carrier transport in 2D bi- 
layer graphene (BLG) has attracted considerable atten- 
tion 0, In BLG, the carriers tunnel quantum me- 
chanically between the two layers leading to a modified 
band dispersion which is approximately parabolic with an 
effective mass of about 0.033me 0. BLG transport thus 
involves dynamics of chiral, parabolic dispersion carriers 
in the zero band gap situation in contrast to the chiral, 
linear dispersion Dirac carrier system for SLG. In ad- 
dition to the considerable fundamental interest in BLG 
transport, there is a great deal of technological interest in 
understanding BLG transport because of the possibility 
of opening a band gap in BLG via an external electric 
field The existence of such a band gap coupled with 
the ability to tune the carrier density using an external 
gate voltage leads to the possibility of graphene-based 
transistor devices with obvious implications for micro- 
electronic and nano electronic applications. 

In this Letter we develop the first comprehensive and 
quantitative theory for BLG carrier transport in the 
presence of random charged impurity (i.e^ long-range 
Coulomb [H]) and short-range scattering Q. Coulomb 
scattering has emerged as the leading, if not the dom- 
inant, scattering mechanism controlling SLG carrier 
transport. In particular, the observed linear dependence 



of the SLG conductivity on carrier density away from the 
charge neutral Dirac point as well as the observed weak 
temperature dependence of conductivity are both nat- 
urally explained as arising from the screened Coulomb 
scattering of chiral Dirac carriers by random charged 
impurities invariably present in the graphene environ- 
ment [7,]. Although the precise quantitative importance 
of charged impurity scattering in determining all aspects 
of SLG transport properties is still occasionally debated 
Q in the literature, there is wide consensus that SLG 
transport both around and away from the Dirac point is 
determined mostly, if not entirely, by screened Coulomb 
scattering from random charged impurity centers. It is, 
therefore, reasonable to assume that charged impurity 
scattering also dominates BLG transport, and to investi- 
gate theoretically the effect of charged impurity scatter- 
ing on BLG conductivity. 

We start by providing a qualitative conceptual discus- 
sion of BLG transport properties vis a vis SLG transport. 
Because of the linear (SLG) versus quadratic (BLG) en- 
ergy dispersion in the two cases, the density of states 
is linear (constant) in SLG (BLG), leading to very dis- 
tinct screening behavior in two cases 9, 10]: q^^p^ (x kp', 
=constant, where qxF is the long wavelength 
Thomas- Fermi screening wave vector and kp oc \[n is 
the Fermi wave vector (with n as the tunable 2D car- 
rier density). Therefore, the screened Coulomb impurity 
potential, Ui{q), behaves very differently in the two sys- 
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{qTF + kp)~^, where 
qxF = qTF^ — i'^Tne^) / {nh^) w lO'^cm""'^, where k is 
the background environmental dielectric constant. Thus, 
the nature of screened Coulomb scattering is qualita- 
tively different in the two systems — in fact, the screened 
Coulomb disorder in the SLG behaves as unscreened 
Coulomb interaction (i.e., v{q) ~ \/q) [§] whereas the 
screened Coulomb disorder in the BLG behaves similar to 
the 2D screened Coulomb interaction although there are 
some differences [lo| . It is, therefore, quite puzzling that 
experimentally the two systems have very similar carrier 
transport properties 0, 13| . Away from the charge neu- 



tral point, CNP, (n = 0), both manifest a conductivity, a, 
approximately proportional to n, whereas for n ~ the 
conductivity is approximately constant (ct ~ 1 — 10 /h), 
forming the much-discussed graphene minimum conduc- 
tivity plateau in both cases. The striking similarity of the 
experimentally observed density dependent conductivity 
a{n) in the two systems presents a serious challenge to 
transport theories since the nature of screened Coulomb 
disorder in the two systems is fundamentally different. 

In this Letter we propose that BLG carrier transport 
is controlled almost equally by two distinct and indepen- 
dent physical scattering mechanisms: Screened Coulomb 
disorder due to random charged impurities in the envi- 
ronment and a short-range disorder (of unknown origin) 
which is important throughout (i.e. all the way from the 
CNP manifesting the plateau-like minimum conductiv- 
ity behavior to the high-density conductivity increasing 
linearly with carrier density). We find that neither pure 
Coulomb disorder nor pure short-range disorder by it- 
self can explain both the density and the T dependence 
of the experimentally observed BLG conductivity behav- 
ior. Earlier theoretical work on BLG transport [H*], con- 
sidering only screened Coulomb disorder, included dras- 
tic approximations (e.g. complete screening) which are 
both unreliable and uncontrolled. We emphasize that 
the actual amount of Coulomb and short-range disorder 
we need to qualitatively and semi-quantitatively explain 
the existing BLG transport data is comparable to that 
already used in understanding the SLG transport prop- 
erties, and as such, there is no arbitrary data fitting pro- 
cedure in our theory. The new element of physics, not 
appreciated in the earlier literature, is that BLG screen- 
ing is quantitatively much stronger than SLG screening, 
rendering the effect of Coulomb scattering relatively less 
important in BLG (compared with SLG). This then leads 
to any other scattering mechanism present in the system 
to become much more important for BLG than for SLG, 
and hence short-range scattering is much more significant 
for BLG than for SLG. This is easily seen by noting that 
the ratio V of screened Coulomb disorder in SLG to BLG 
goes as r = uf^'^/uf^'^ = (I + ^tf/^f) « {1 + 11.5/ yM) 
where n — (n/10^^). Thus, r 3> 1 in the typical graphene 
experimental carrier density regime. 

The density and temperature dependent BLG conduc- 
tivity is given within the Boltzmann transport theory by 
a = (eVm) / de7^(e)er(e) (-df/de) , where D{e), r(e), / 
are respectively the energy dependent density of states, 
the transport relaxation time (which depends explicitly 
on the scattering mechanism), and the (finite-T) Fermi 
distribution function, r is given by 




X[l - COSeikk-lfJ (Csk - Esk') ; (1) 

where esk = sh^k"^ /2ra is the carrier energy for the 
spin/pseudospin state 's' and 2D wave vector k, (y^k.sk') 
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FIG. 1: (Color online) The wave vector dependent 2D BLG 
polarizability 11 (in units of the density of states at the Fermi 
level Nf) for different values of T/Tf =0, 0.1, 0.5, 1.0, 2.0, 
where Tf ~ EF/ks is the Fermi temperature (Tf ~ 831K for 
n = 10^^ cm~^) and the dielectric function is given by e = 
1-1- Dell where Vc is the 2D Coulomb interaction. The ordinate 
can be taken as a measure of the strength of BLG screening 
in units of Thomas- Fermi screening constant qtf ~ lO^cm"^. 



is the matrix element of the appropriate disorder poten- 
tial, g{6) = [1 + cos(26')]/2 is a wave function form factor 
associated with the chiral matrix of graphene (and is de- 
termined by the BLG band structure). In Eq. ([1]), no is 
the appropriate 2D areal concentration of the impurity 
centers giving rise to the random disorder potential V. 
We consider two different kinds of disorder scattering: {i) 
screened Coulomb disorder Ui, and (ii) short-range disor- 
der Vb, writing no|(V;k,sk')P = "•il^iP = ni\vi{q) / e{q)\'^ 
andno|(^sk,sk')P = ndV^, respectively. Here g = |k-k'|, 
Vi{q) — (27re^/K(7)e^''* is the 2D Coulomb potential, and 
Vq a constant short-range (i.e. a (5-function in real space) 
potential. (Note that we use Ui to denote the charged 
impurity density.) In the expression for Vi, d is the aver- 
age distance of the charged impurities from BLG. We find 
that in contrast to SLG, in BLG the conductivity is essen- 
tially independent of d for d < 10 A due to the short BLG 
screening length. The dielectric screening function e{q) 
entering the effective screened Coulomb disorder, which 
depends on both q and T, was calculated in Ref. [l^] at 
T — 0, and we have now generalized it to T ^ as shown 
in Fig. 1. Note that BLG screening peaks at 2kp which 
is also the important BLG scattering wave vector due to 
the form factor g in Eq. (1). 

It is straightforward to calculate the analytical density 
dependence (for T/Tp <C 1) of BLG conductivity from 
the above formulae: (j{n) ^ for unscreened Coulomb 
disorder; cr(n) ~ n for overscreened Coulomb disorder; 
(j{n) ^ n", 1 < a < 2 for screened Coulomb disorder; 
a{n) ^ n for short range disorder. We emphasize that 
a{n) is a density dependent exponent which varies slowly 
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FIG. 2: (Color online) (a) Density dependence of BLG con- 
ductivity for different temperatures, T=0, 50, 100, 150, 200, 
300 K (from bottom to top), rii = 10^^ cm~'^ and UdVQ — 
2.0(eV^A)^. Top inset shows ct as a function of T in pres- 
ence of short range disorder. The scaled conductivity applies 
for all densities. Bottom inset shows cr as a function of T in 
presence of screened Coulomb disorder for different densities 
n = 5, 10, 30 X 10" cm~^ (from bottom to top). 



changing from 1 at low density to 2 at high density - in 
the BLG experimental density range a « 1.2. Increasing 
temperature, in general, suppresses screening, leading to 
a slight enhancement of the exponent a. 

It is obvious from the above discussion that Coulomb 
disorder by itself cannot explain the experimentally ob- 
served linear density dependence of cr(n) in BLG. By con- 
trast, the SLG conductivity for Coulomb disorder, either 
screened or unscreened, follows the linear a{n) ^ n be- 
havior whereas the short range disorder leads to a density 
independent cr(n) n° behavior. Thus, the carrier trans- 
port physics is fundamentally and qualitatively different 
in SLG and BLG: SLG is always dominated by Coulomb 
disorder except at extreme high densities, whereas in 
BLG the short-range disorder is always equally, if not 
more, important for all densities. 

In Fig. [2] we present our full numerically calculated 
finite T BLG conductivity as a function of n for rea- 
sonable representative values of disorder strength. The 
results shown in Fig. [2] are in good agreement with recent 
experimental results [3, Q, but we do not make any at- 
tempt to obtain exact quantitative agreement since the 
experimental results show substantial sample-to-sample 
variations Instead we discuss the salient qualitative 
features of our results: (?) The calculated density depen- 



dence is approximately linear over a wide density range 
as seen experimentally; (ii) the temperature dependence 
is very weak at higher densities as observed in recent ex- 
periments 0; [iii) at low densities, where T/Tp is not 
too small, there is a strong insulating-type T dependence 
arising from the thermal excitation of carriers (which is 
exponentially suppressed at higher densities) and energy 
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FIG. 3: (Color online), (a) n(r) of BLG at the CNP for a 
single disorder realization with rti = lO^^cm"^ and d = 1 nm. 
(b) Disorder averaged P(n), at the CNP for BLG (SLG) red 
(blue) for m = lO^cm'^ and d = 1 nm. For SLG P(n = 0) ^ 



0.1, out of scale. The corresponding rir, 
BLG and 1.2 lO^cm-^ for SLG. 



is 5.5 lO^cm"^ for 



averaging, as observed experimentally 0]; (w) when the 
dimensionless temperature is very small {T /Tp <^ 1) our 
theory necessarily predicts (see lower inset of Fig. [5]) a 
weak linear-in T metallic T dependence arising from the 
temperature dependence of the screened charges impurity 
scattering, i.e. the thermal suppression of the 2fci?-peak 
associated with back-scattering in Fig. [T] ~- this effect is 
enhanced in BLG due to the importance of backscatter- 
ing whereas it is suppressed in SLG. By contrast, for the 
short-range disorder a always increases with T, as shown 
in the upper inset of Fig. [51 

While the calculations for Fig. [5] are all carried out for 
graphene on Si02 substrate (corresponding to k w 2.5), 
we have carried out detailed calculations of cr(n, T, k) as 
a function of the effective background dielectric constant 
(not shown) since screened Coulomb disorder should 
manifest a strong dependence on k. Our calculations 
show the K-dependence of the net conductivity to be al- 
most unobservable - for example, covering BLG on Si02 
with ice, thus changing k from 2.5 to 3.5, would only in- 
crease the conductivity by 1%, for rii — lO^^cm"^, and 
10%, for = 10i2cm-2. 

The Boltzmann theory presented so far assumes a com- 
pletely homogeneous carrier density landscape over the 
BLG sample, and leads to a{T = 0) at the CNP 
(n = 0). It is, however, well known that SLG breaks up 
into an inhomogeneous landscape of electron-hole pud- 
dles around the CNP [H, H [III ■ We expect the same 
physics of electron-hole puddles to dominate the BLG 
properties around the GNP, and this will give rise to a 
finite "minimum conductivity", CTmin = cr(n = 0) even in 
the limit T ^ 0. We have investigated the properties of 
BLG puddle formation by calculating the ground state of 
the BLG in the presence of random charged impurities 
solving the density functional equations within a local 
density approximation. [l3| Using the resultant ground 
state density landscape n(r) with n = {n{r))r and aver- 
aging over disorder realizations we have then calculated 
the BLG transport properties using the effective medium 
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FIG. 4: (Color online). Conductivity as a function of 
n obtained using the EMT for Ui — lO^^cm"^, UdVo ~ 
2 (eVA)^ and several values of T, from top to bottom: 
T = 300, 200, 100, 0. From the slope for n > lO^^cm"^ 
for r = we extract a mobility of 4000 cm^/Vs. In the lower 
(upper) inset (T,nin (cin ~ 1-9 lO^'^cm"'^)) as a function of T 
for the same values of disorder strengths. 



theory, EMT [15| . The details of the calculations will be 
presented elsewhere, here we present the realistic conduc- 
tivity results in the puddle dominated regime obtained 
using our EMT. 

In Fig. [3] (a) we show our calculated zero-density pud- 
dle structure for the BLG for a single disorder realization 
with Tii = lO^^cm^^. We note that by neglecting the dis- 
order induced inhomogeneities BLG (SLG) would have 
perfect (vanishing) linear screening properties as n — > 0, 
whereas the strong carrier density fluctuations and as- 
sociated electron-hole puddle structure induced by the 
charged impurities close to the CNP are qualitatively 
very similar in BLG and SLG. The linear (SLG) versus 
parabolic (BLG) carrier dispersion has the quantitative 
effect of modifying the form of the probability distribu- 
tion P{n) with the result that P{n) is much wider in 
BLG than in SLG, as shown in Fig. [3] (b). 

Finally in Fig. [4] we show cr{n) for BLG for several 
values of T, taking into account the inhomogeneity of 
the puddles. These results are calculated numerically 
through the EMT by combining the density functional 
electronic structure calculation with the full Boltzmann 
transport theory. At high density (3> n-i), the theory 
gives the same results as that obtained from the Boltz- 
mann theory in the homogeneous case (i.e. Fig.[2|), but at 
low densities there are significant deviations. For exam- 
ple, we now get a well-defined finite amin even at T = 0. 
We show the T dependence of the conductivity for n = 
and rii in the insets obtaining good agreement with 
the experimental finding of a strong insulating T depen- 
dence of iTniin- Our Calculated CTmin depends weakly on Ui 
with no universally discernible functional dependence on 
Ui. For the range of values of n and rii considered using 



a more realistic hyperbolic BLG band dispersion in the 
theory does not change our results and conclusions in any 
qualitative manner. A small BLG band gap (^ Ep) at 
k — also has no significant effect. 

We conclude by emphasizing the similarity and the 
difference between BLG and SLG transport from the 
perspective of our transport theory considerations. We 
find that both manifest a non-universal sample depen- 
dent minimum conductivity at the CNP arising from the 
electron-hole puddle formation due to the presence of 
long-range Coulomb disorder. We find that BLG mani- 
fests a fairly strong insulating temperature dependence of 
Cmin, which is consistent with experimental observations. 
We find that at high density both BLG and SLG manifest 
a linearly increasing conductivity with increasing carrier 
density, as observed experimentally. However the physi- 
cal origin for this linear dependence is quite different in 
the two systems: while in the SLG this linearity arises 
entirely from Coulomb disorder, in the BLG we must in- 
voke the existence of a short-range disorder scattering 
to explain the linearity. The importance of short-range 
scattering in BLG compared with SLG is understand- 
able based on BLG screening being much stronger than 
SLG screening. The existence of the short-range disor- 
der naturally explains the strong insulating temperature 
dependence of BLG conductivity since short-range disor- 
der leads to a cr(T) increasing with T for all T, whereas 
Coulomb scattering always leads to a (t{T) decreasing 
(linearly at first) with increasing T for T <^ Tp. A pos- 
sible direct way of experimentally investigating the com- 
peting effects of short-range versus Coulomb scattering 
in BLG will be to plot the measured low-temperature 
conductivity against the carrier density on a log-log plot 
to see whether the power law exponent at higher den- 
sities is larger than unity, thus implying the dominance 
of Coulomb scattering. The origin of the strong short- 
range disorder in BLG has to be considered unknown 
at this stage. It could be neutral defects such as lat- 
tice vacancies (only a minuscule concentration of 10^^ or 
less would be needed), or ripples, or resonant short-range 
scattering, or some other mechanism. 
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